Numerical Simulation of Detonation in Condensed Phase Explosives

This report describes the development of a two-dimensional multi-material Eulerian hydrocode to model the effects of detonating condensed phase explosives on surrounding materials. The code solves the Euler equations for the conservation of mass, momentum, and energy for an inviscid, compressible fluid. Operator splitting is used to reduce the two-dimensional calculation into coupled one-dimensional equations, which are then solved using the Flux-Corrected Transport (FCT) algorithm of Boris and Book. Non-reacting materials are described using either a perfect gas, MieGruneisen, or Tait equation of state, while the energetic materials are described using either a BKW equation of state and Forest Fire reaction rate model, or the JWL equation of state and the Ignition and Growth reaction rate model. A modified Young’s algorithm is used to maintain a sharp interface between different materials on the computational mesh. A brief description of the major components of the coding is provided, and then several applications of the code are described, including the simulation of bullet impact experiments, the underwater sympathetic detonation test, and the modified gap test.